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ABSTRACT 


Diffraction  theory  is  used  in  analyzing  the  radiation  characteristics 
of  typical  horn  antennas.  The  far- side-lobe  and  bacK-lobe  radiation  has 
been  solved  without  employing  equivalence  principles  which  are  impractical 
in  the  problem. 

A  comer  reflector  with  a  magnetic  line  source  located  at  the  vertex 
is  proposed  as  a  model  for  the  E-plane  radiation  of  horn  antennas.  Dif¬ 
fraction  theory  is  applied  to  formulate  the  radiation  pattern  of  this  antenna 
model.  A  complete  pattern,  including  multiple  interactions  and  images 
of  induced  line  sources,  is  obtained  in  infinite  series  form.  Diffraction 
mechanisms  are  used  for  appropriate  approximations  in  the  computations. 
Emphasis  is  made  on  the  continuity  of  the  radiation  patterns.  The  com¬ 
puted  patterns  are  in  excellent  agreement  with  measured  patterns  of 
typical  horn  antennas.  Radiation  intensity  of  the  back-lobe  relative  to 
main- lobe  intensity  is  obtained  as  back-to-front  ratio  and  plotted  as  a 
function  of  antenna  dimensions.  A  set  of  measured  back-to-front  ratios 
is  presented  to  demonstrate  the  validity  of  computation. 
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THE  E- PLANE  RADIATION  PATTERN  OF  AN 
ANTENNA  MODEL  FOR  HORN  ANTENNAS 


I.  INTRODUCTION 

Electromagnetic  horns  are  widely  used  either  as  feeds  or  direct 
radiators  in  antenna  systems.  Therefore,  a  complete  and  practical 
understanding  of  the  radiation  mechanisms  is  important.  The  purpose 
here  is  to  provide  this  understanding  by  using  diffraction  theory  as  a 
tool  fcr  analyzing  radiation  characteristics  of  horn  antennas. 

Experimentally  radiation  patterns  of  different  horn  lengths  and 
flare  angles  in  both  E-  and  H-planes  have  been  measured  by  Rhodes[l], 
These  patterns  are  adequate  only  in  the  vicinity  of  the  main  lobe.  For 
the  regions  of  the  far  side-lobes  and  back  lobes,  which  are  important 
in  EMC  (Electromagnetic  Compatibility)  problems,  more  accurate 
measurements  are  required, 

Analytically,  the  descriptions  c£  the  propagating  modes  in  a  horn 
are  summarized  by  Kraus  £  2].  In  general,  aperture  techniques  must 
be  used  to  calculate  the  radiation  pattern,  and  it  is  assumed  that  the 
aperture  distribution  is  that  of  the  incident  wave  and  is  zero  outside 
the  aperture.  The  patterns  thus  obtained  give  satisfactory  main  lobes 
and  near  side-lobes.  As  for  far  side-iobes  and  back  lobes,  Schelkonof^s 
equivalence  principles  cculd  be  applied  if  current  distributions  on  the 
outer  surfaces  were  known.  For  a  horn  antenna,  extreme  difficulty  is 
involved  in  accurately  describing  the  current  distributions.  Furthermore, 
difficulty  would  arise  in  evaluating  the  consequent  surface  integrals.  The 
impracticality  of  the  equivalence  principles  leaves  radiation  problems  of 
the  far  side-lobes  and  back  lobes  still  unsolved. 

Because  of  the  increasing  power  levels  of  individual  antennas  and 
the  increasing  numbers  of  radiating  systems  in  particular  areas,  the 
characteristics  of  far  side-lobes  and  back  lobes  have  become  increasingly 
important  in  EMC  problems.  Since  these  radiation  characteristics  are 
caused  by  diffraction  of  electromagnetic  waves,  it  has  been  a  natural 
consequence  in  recent  years  to  apply  diffraction  theory  to  the  radiation 
problems  of  horn  antennas. 
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In  1962,  Kinber  [3]  derived  horn  patterns  and  the  coefficient  of 
coupling  between  two  adjacent  horns  by  diffraction  theory,  Examples  are 
given  for  both  E-plane  and  H-plane  patterns,  in  which  discontinuities  are 
pointed  cut,  and  emphasis  has  not  been  made  on  side  and  back  lobes.  In 
1963  Ohba[4j  used  diffraction  theory  to  compute  the  radiation  pattern 
in  the  H-plane  of  a  corner  reflector.  Disagreement  was  noted  as  a  re¬ 
sult  of  neglecting  contributions  from  the  other  edges  of  the  corner  re¬ 
flector,  In  1964,  Russo,  Rudduck,  and  Peters[5)  employed  diffraction 
theory  with  proper  assumptions  to  obtain  E-plane  patterns  of  a  thin- 
edged  and  a  thick-edged  horn.  Only  the  first-order  diffraction  terms 
were  used  to  compute  the  thin-edged  horn  patterns.  The  results,  with 
possible  discontinuities  left  in  the  side  and  back  lobes,  are  in  good 
agreement  with  the  measured  pattern.  Even  though  the  higher-order 
diffraction  farms  and  the  reflections  inside  the  horn  are  neglected,  the 
combination  of  the  employed  concepts  and  assumptions  constitutes  a 
new  method  of  analysis  for  horn  antennas.  We  shall  follow  this  new 
method  to  develop  a  more  complete  analysis  by  including  the  previously 
neglected  higher-order  diffractions  and  the  reflections  inside  the  horn 
antenna. 


II.  RADIATION  MECHANISMS 

The  proposed  antenna  model  ia  a  corner  reflector  formed  by  two 
perfectly  conducting  plane  walls  intersecting  at  an  angle,  26g,  as  shown 
in  Fig.  1.  If  we  let  the  comer  reflector  be  infinite  in  extent  along  the 
s-direction,  the  problem  is  thus  reduced  to  a  two-dimensional  one.  The 
primary  source  is  a  magnetic  line  current  assumed  at  the  vertex  S. 

This  assumption  considerably  simplifies  treatment  in  the  principal  E- 
plane  of  a  horn  antenna  fed  by  a  waveguide  supporting  the  TE|*  mode. 

The  angular  coordinate,  %  shown  in  Fig.  1,  is  the  common  refer¬ 
ence  angle.  Tn-t  angles  *{«  are  the  diffraction  angles  referred  to  each 
individual  wedge  at  which  diffraction  occurs.  The  angles  a  are  called 
incident  angles  of  illumination.  All  the  first  subscripts  refer  to  the 
points  at  which  diffraction  occurs,  while  the  second  subscripts  refer 
to  the  points  of  origin  of  incident  rays.  This  notation  will  be  used  through¬ 
out  the  following  discussions.  There  are  four  wedges  (A,  B,  S,  and  W)  to 
be  treated  by  diffraction  theory.  Wedges  A  and  B  have  sere  wedge  angle, 
while  wedges  S  and  W  have  wedge  angles  of  and  2(w*6E),  respectively. 
The  property  of  symmetry  of  the  reflector  will  be  used  to  simplify  ths 
problem  by  considering  only  the  upper  half  of  the  pattern. 

In  diffraction  theory,  a  uniform  cylindrical  wave  is  radiatsd  from 
the  primary  source  in  the  region  -0g<  8<  This  uniform  cylindrical 
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wave  is  called  the  geometrical  optics  wave  which  illuminates  wedges 
A  and  B.  The  diffractions  at  A  and  B  caused  fay  this  illumination  are 
called  the  first-order  diffractions  which  have  the  solutions  shown  in 
Eq.  (32)  (Appendix  A),  The  solutions  are  directional  cylindrical  waves 
radiated  from  the  wedges.  The  geometrical  optics  rays  from  S  and  the 
first-order  diffracted  rays  from  A  and  B  are  shown  in  Fig.  2a.  The 
first-order  radiation  pattern  in  the  far-field  can  now  be  obtained  by 
superposition  of  the  far-field  intensities  of  the  primary  source  and  the 
two  induced  sources. 

To  consider  the  diffraction  process  further,  one  can  observe  from 
Fig.  2b  that  the  induced  source  at  B  illuminates  wedges  A,  S,  and  W  to 
give  three  second-order  diffraction  terms.  In  the  same  manner,  the 
first-order-induced  source  at  A  illuminates  wedges  B,  S,  and  W  to  give 
three  more  second-order  diffracted  waves.  These  six  second-order- 
induced  sources  will  continue  to  give  third-  and  higher-order  diffraction. 
The  induced  intensity  becomes  smaller  with  increasing  order,  and  the 
phase  delay  of  successive  illumination  can  be  properly  taken  into  account. 
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Fig,  2.  Radiation  mechanism  of  the  antenna  model,  a)  Direct  rays 
and  the  lst-order  diffracted  rays  due  to  illumination  from 
primary  source  at  S;  b)  The  2nd-order  diff  .  actions  due  to 
the  lst-order  illumination  from  B;  c)  The  first  image  in  the 
lower  wall  due  to  the  lst-order  diffraction  at  A;  and  d)  The 
2nd-order  diffraction  due  to  the  first  image  in  the  lower  wall. 

Since  the  reflector  has  perfectly  conducting  walls,  the  diffracted 
waves  from  A  and  B  are  reflected  by  the  walls.  The  first-order  diffracted 
rays  from  A  are  partially  reflected  by  the  "lower  wall,  as  shown  in  Fig. 

2c.  The  reflected  rays  can  be  described  by  the  image  method  of  Appendix 
B.  Figure  2d  shows  that  wedge  A  is  illuminated  by  one  of  the  first-order 
images  from  the  lower  wall.  The  number  of  images  is  determined  by  the 
flare  angle  of  the  reflector.  The  effects  of  the  reflector  wadis  can  then  be 
taken  into  account  by  the  images  and  the  subsequent  diffraction  of  the  images. 

When  the  process  of  diffraction  and  reflection  described  above  is 
completed,  the  far-field  patterns  of  th '  reflector  antenna  can  be  obtained 
by  superimposing  the  contributions  from  the  primary  source  at  S;  the 
induced  sources  at  A,  B,S,  and  W;  and  the  images  in  both  falls.  Form¬ 
ulation  of  the  pattern,  including  all  orders  of  diffraction  and  reflection, 
is  obtained  in  the  next  section. 
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III.  FORMULATION  OF  SOLUTION  IN 
INFINITE  SERIES  FORM 

Diffraction  theory  is  briefly  treated  in  Appendix  A  for  both  plane- 
wave  and  uniform  cylindrical-wave  illuminations.  The  wedges  of  the 
corner  reflector  studied  here  are  all  illuminated  by  cylindrical  waves. 

In  Appendix  B,  a  method  of  images  is  presented  to  describe  the  image- 
waves  formed  inside  the  antenna.  The  diffraction  theory  and  the  method 
of  images  are  used  in  the  following  to  formulate  a  solution  for  our  problem. 

First,  referring  to  Figs.  1  and  2a,  the  geometrical  optics  wave 
radiated  from  the  primary  source  S  is  a  uniform  cylindrical-wave, 
defined  as 


(1) 


v*(0)  =  1 


-0e<  0<  +  8e* 


where  v*  has  point  S  as  phase-reference,  and  outside  the  defined  region 
v*  is  identically  zero.  The  cylindrical-wave  propagation  factor  R“2 
Exp.  (-jkR)  to  the  far-field  is  suppressed  in  Eq.  (11)  because  only  the 
angular  dependence  is  of  interest. 

Wedges  A  and  B  are  illuminated  by  the  cylindrical  wave  from  S 
with  zero  incident- angle.  Since  there  is  no  reflection  term,  tne  dif¬ 
fracted  waves  from  A  and  B  have  only  one  term  of  the  solutions  given 
by  Eq.  (32).  Excluding  the  portion  ci  waves  diffracted  into  the  corner 
reflector,  the  waves  directly  diffracted  to  the  far-field  can  be  written 
as 

(t)  _ 

(2)  ^as *  vb(Pjt»  w  ~0e  ^  0*  ^)  ~  y  — iS  n  * 


Dgg  =  vb(pe,  w  -  0£  -  6,  2)  J  —  +  0e)  » 

^AS  =  ‘l' AS  =  *-0E  +  9' 

^BS  s  ^BS  =  *~SE  -  0. 
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(2) 

(cont) 


SA  = 

tJi:  v&/|  <4*4Si 


nA  =  nB  =  2» 


where  and  DBg  designate  diffraction  at  A  and  E  because  of  illumi¬ 
nation  from  S.  The  superscript  {1)  means  first-order  diffraction.  The 
expressions  of  <^s  in  terms  of  6  can  be  obtained  from  Fig.  1.  The 
argument  n  is  equal  to  2  for  both  A  and  B  because  they  have  zero  wedge 
angle.  It  should  bt  noted  that  the  notation  vg  follows  the  form  of  the 
original  solution  and  the  subscript  B  has  no  connection  with  the  wedge  B. 

The  first-order  radiation  patterns,  neglecting  the  reflections  in¬ 
side  the  comer  reflector,  can  now  be  obtained  by  simply  superimposing 
the  terms  in  Eqs.  <1 )  and  (2),  The  discontinuities  in  V  at  6=  +  Sgin 
Eq.  (1)  are  eliminated  by  and  respectively,  in  the  manner 

illustrated  in  Eq,  (34).  Although  the  pattern  is  made  continuous  at 
6  =  ±0£,  two  sets  of  new  discontinuities  at  8  =  ±rr/2  and  i(Tr  +  8g)  are 
observed  in  Eq.  (2).  Therefore,  the  first-order  pattern  in  general  has 
discontinuities  at  these  directions. 


Let  us  next  examine  the  reflections  of  the  first-order  rays  dif¬ 
fracted  into  the  reflector.  Since  the  diffracted  waves  from  A  and  B  are 
symmetrical  with  respect  to  8  =  0,  as  can  be  seen  from  Eq.  (2),  the 
image-waves  from  the  reflector  can  be  treated  as  in  Appendix  B.  The 
images  are  formed  symmetrically  in  the  lower  and  upper  walls.  Re¬ 
placing  U  of  and  in  Eq.  (2)  by  (-2i&g  -  8)  and  (2i8£-  8),  respec¬ 

tively,  the  image-waves  from  two  walls  can  De  obtained  similar  to  Eq. 
(44)  as 


(3) 


(IL^i  =  vB{pE,7r-(2i+l)eE-8,2),  f  ~(i+1)8E£eif  ~i0E  : 


(I 


(i) 

U)i  =  vB(PE.  w  -(2i+1)9e  +  2)» 


I-(Hl)9J>9>  - 


T-% 


i  =  1,  2,  3 


-h;  and 


h  (the  largest  integer)<  -J-, 

-  2% 
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where  the  subscripts  L  and  U  indicate  that  the  image  terms  are 
the  lower  and  the  upper  wails,  respectively.  The  number  of  images  in 
each  wall  is  equal  to  h,  as  discussed  in  Appendix  B.  The  division  of 
regions  is  also  studied  in  Appendix  B.  When  the  ratio  ir/ZO^is  not  an 
integer,  the  valid  region  for  the  last  images  should  be  modified  to 


{4) 

and 


^•-{h+l)Sg<  0<  w-{2h+l)0^  for  the  lower-wall, 


—•  -(h+l}8£ 


>  0>  -  fir  - (Zh-t-1  for  the  upper- wall  . 


Each  term  in  Eq.  (3)  has  its  properly  defined  regions  and  is  set  zero 
outside  the  region.  The  first  image  in  the  lower  wail,  caused  by  the 
1'irst-order-diffracted  rays  from  A,  is  shown  in  Fig,  2c.  Figure  3  shows 
the  images  in  the  lower  wall  for  h  =  4.  It  is  noted  that  the  true  image 
waves  of  the  diffracted  waves  from  A  are  those  with  i  odd  in  the  lower 
wall  and  i  even  in  the  upper  wall. 

Descriptions  of  the  first-order  diffracted  waves  from  A  and  B  and 
their  reflected  waves  have  been  completed  above.  The  higher-order 
terms  to  be  treated  in  the  following  discussion  are  necessary  for  cases 
in  which  small  dimensions  are  encountered  or  high  accuracy  is  desired. 
Physically,  the  higher-order  terms  describe  the  effects  of  illumination 
of  edges  by  the  lower-order-induced  sources  and  their  images.  Mathe¬ 
matically,  they  are  required  to  overcome  the  discontinuities  of  the 
lower-order  terms  in  the  radiation  pattern.  Talking  the  two  first-order 
diffraction  terms  in  Eq.  {2},  for  instance,  the  discontinuities  mentioned 
earlier  can  only  be  eliminated  by  taking  into  account  the  second-order 
diffraction  •'  i  the  specified  directions. 

At  9  =  ir/2  in  Eq.  (2),  w-edge  A  is  illuminated  by  the  first-order- 
induced  source  at  B.  This  intensity  of  illumination  from  B  to  A  is 
called  the  first-order  coupling  coefficient. 


as  shown  in  Fig.  2b.  Because  of  symmetry,  the  first-order  coupling 
coefficient  from  A  to  B  in  the  direction  9  =  -n/2  is  equal  to  C^jg.  There¬ 
fore,  using  Eq„  (2),  we  have 
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Fig.  3.  Four  image*  in  the  lower-wail. 


'AB  ~  CBA  =  vB(f>E'  §  '  9e*  2). 


re.pectively^^J^ :°^plingl^efffcient«>  cai^ bt  **  9  =  *(*  +  9e)* 
manner  a*  *  «i.cient*  can  be  obtained  in  the  .ame 


c*  •  -  r~ 

^ur  a  -  C, 


'B^E*  2*»  2) 


Since  the  diffracted  wave*  , 

neighborhood  of  a  certain  angle,  it  i.  a  Zo £VaryinS  functions  in  the 
A,  B,  and  W  are  illuminated  by  uniform  fvif  ffr°xima^°n  that  wedge, 
sitie.  shown  in  Eq..  (5)  and  {T>.  UnZ  til  ^  W*V”  °f  the  Re¬ 
order  diffracted  wave,  can  be  obtained  a  1  a8*UniPtioa’  *e  .econd- 
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=  C^B  [vB(b.  |w,  2)+  vB!b,  *1  -29e+8,  2)j, 

"  f  i  91  *  +  ®E  * 

DBA  ■  CfiA  -6-  2>  +  -2«E-8.  2)1  . 


|  >  e>  -  {*  +  %)  ; 


DWAa  CWA  QE±Q*2w~  °E‘>  and 


d$b  =  C^B  [vg{p£,  2 ir  -  0E-  6,  nw)],  -6£  >  8>  -{2»  -  8g)  , 

where  the  arguments  used  can  be  obtained  from  Fig.  1  as 


*A£ 

=  di  +  =  (* 

yAB  AB 

“V 

♦* 

BA 

=  iL  ^  ®  ,  =  (tt 

^BA  BA 

“  ®E 

*WA 

=  4<wa  =  6  -  eE» 

d> 

WB 

=  +WB  '  21  •  s£  - 

6, 

nw  = 

26-, 

2  ”  — ^  . 

rt 

I“6e>’ 


where  nyy  is  obtained  from  wedge  angle,  (E-n^y)*,  of  wedge  W.  Next, 
we  note  in  Eq.  (3)  that  the  image  waves  are  discontinuous  at  angles 


6=  t  (j  -  %>•  1  *  1.2.3  -  -  (h-1)  , 
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at  which  the  wedges  A  and  B  are  illuminated  by  the  rays  from  image? 
of  lower  and  upper  walls,  respectively.  Figure  3  shows  the  geometry 
of  the  images  with  h=4  in  the  lower  wall.  The  coupling  coefficients 
from  the  images  to  wedge  A  can  be  obtained  from  Eqs.  (3)  as 


CAi  =  (e  =  §  -  i8E  >*  1  =  l»  Z> . ^h-1>* 

By  symmetry,  the  coupling  coefficients  from  the  images  in  the  upper 
wall  to  wedge  B  are  equal  to  C^({  as 


(9) 


=  CB.  -  vb^e»  "T  “  2),  i=l,  2,  -  -  (h-1). 


The  second-order  diffraction  at  A  and  B  illuminated  b;  the  images  can 
now  be  written  as 


(10) 


DS  =  Ca!  [VB(  Pi-  5  +  i0E+e-  2>  +  VB(  Pf ' T -<i+2>V0-  2l]  - 


and 


_  r(a) 

Bi  -  CBi 


-|<e<a+eE, 


rB(  pi‘  2  +  lV8*  2>  +  VPl'T-< i+2)0E*e,  2)j 


|>  0>  -  (»  +  9e)p 


i  =  1,2, - (h-1), 


where  tne  arguments  can  be  obtained  by  using  Fig.  3  with  h  =  4  as  reference. 
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(11)  Pi  =  Pi-i  COS  0E  T  pD  cos  iO£j  and 

(cont) 

i  =  1,2, - (h-1). 


The  second-order  diffraction  terms  obtained  in  £q.  (iO)are  appropriately 
arranged  for  each  boundary  of  the  defined  regions  in  Eq.  (3),  except 
that  the  last  boundary  is  given  by  8  =  +  (tt/2  -(h+l)©^)  if  ir/2©g  is  an  inte¬ 
ger,  or  6  =  ±  (jr-(2h+i)0E)  if  ir/26 £  is  not  an  integer.  The  boundaries 
in  either  case  correspond  to  the  directions  in  which  wedge  S  is  illuminated 
by  the  induced  sources  at  A  and  B.  Therefore,  the  coupling  coefficients 
from  A  and  B  to  S  can  be  obtained  from  Eq.  (3)  by  symmetry  as 

*-fjA  =  CSB  =  vB^pE*  °’2^  * 

which  gives  rise  to  the  second-order  diffraction  at  wedge  S  as 

t13>  DSA=41i[vB(,’E.Ve-Vi 

DSB  =CSB  9E+9'  3sU 

®SA  =  eE  -  e* 

^SB  =  eE+  9»  and 

nss2-nw  =  -~  • 


Now,  we  have  completed  the  descriptions  of  all  second-order 
diffractions  which  physically  take  into  account  the  effects  of  illumination 
by  the  first-order-induced  sources  and  mathematically  eliminate  all  the 
first-order  discontinuities.  Summing  up  the  second-order  diffraction 
at  A,  B,  W,  and  S  gives 


-f  i9i*+ 


it,  i..w'MIi|!!||I|,,,i'ii;I!",iii  1 !  ®  u  f'tt 


(14) 

cont. 


a  DL2l  + 


h-i 


Y  *>, . 


i~i 


—  >  0  >  -  (tr-f  Or)  : 
L  ~  “ 


DWA 


Ts{l) 

^V/B  * 


Ds  ‘  *  °s  A  +  4t  -  1  el  8e  . 

where  Eqs«  (7),  (10),  and  (13)  can  be  used  for  computation.  Following 
the  same  procedure  used  to  obtain  Eq.  (3),  the  second- order  image 
waves  from  the  lower  and  upper  walls  can  be  obtained  as 


(15) 


vB(b,  f -2i8E‘6*  2)  +  V-b’  ~T  -2<i+l)0E-9,  2)J 

+  I  CAk[vB<P>k>  2>  +  ^-(3»2)9E-e,  2)| 

k=l 


I-(i+l)0E<  Q<1~  %  ; 


(ffii  =  vB(b,  f -2iGE+0.  2)  +  vB(b,-^  -2(i+l)eE+8,  2) 


h-1 

+ 1  CBk 


vB(Pk'  y  +  vB^k*  ”~(3i+2)8j£+0,  2) 


3ir 


f-(i+J)eE]  >6>-[!-<8e 


and 


i  =  i,  2,  3  —  -  h. 
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where  the  8  of  Da  and  d(^  in  Eq.  (14)  have  been  replaced  by  {i^iOjr-Q}, 
respectively.  Note  that  the  boundary  of  the  last  image  i=h  should  follow 
Eq.  (*),  if  tr/ 29j£  is  not  exactly  an  integer. 

We  have  observed  above  that  while  the  first-order  discontinuities 
are  eliminated  by  the  inclusion  of  second-order  diffractions,  new  dis¬ 
continuities  occur  again  at  the  boundaries  of  the  regions  defined  in  Eqs. 

(14)  ana  (15).  These  second-order  discontinuities  can  be  eliminated  only 
by  introducing  third-order  diffraction.  The  higher  the  order  of  diffrac¬ 
tion  the  smaller  will  be  the  magnitude  of  discontinuities.  It  is  theoretically 
possible  to  consider  the  order  of  diffraction  as  high  as  desired.  In 
other  words,  magnitude  of  discontinuities  can  be  made  negligibly  small 
if  sufficient  order  of  diffraction  is  included. 

For  completeness,  let  us  consider  higher-order  diffraction.  The 
third-order  diffraction  at  wedge  A  is  caused  not  only  by  illuminations  from 
the  second-order-induced  sources  at  B  and  the  images  in  the  lower  wall, 
but  also  by  illumination  from  wedges  W  and  S.  In  the  same  manner  as 
Eq.  (14),  the  resultant  third-order  diffraction  at  A  can  be  obtained  as 


DW=  d<3> 
A  AB 


d(’) 

AW 


D(j} 

AS 


Similarly,  by  including  all  the  (m-1  ^-order-induced  sources  and  images, 
the  mth-order  diffraction  at  wedges  A,  B,  W,  and  S  and  the  images  of  the 
corner  reflector  can  be  summarized  as 


<16>  Diml  *  dS!  +  Diw +  DS’ + 1 D 


(m) 

Ai 


i=l 


h-1 


(m)  (m)  (m)  (m)  Vs  (m) 

DB  DBA  +  °BW  +  DBS  +  /  °Bi  * 


BS  / 
i=l 


D(m)  =  D(m)  +  D(m) 
°W  °WA  +  °WB  ' 

(m)  (m)  (m) 

DS  =  +  °SB  ' 
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(16) 

cont. 


(mh  (m) 

(It  h  =  Da  {^-2i6„-9},  and 


(ffii  =  D{Bm)  (6=+2ieE-0)  , 


where  the  valid  regions  are  identical  to  those  defined  in  Eqs.  (14)  and 
(15).  The  components  in  Eq.  (16)  can  be  obtained  analogous  to  Eqs. 

(7),  (10),  and  (13)  as  follows  (the  equal  sign  following  3  in  the  parentheses 
means  "replaced  by"): 

(17)  =  C(™Bl)  [vB(b.f  +9,  2)  + vB(b,  ^-26 e+9,  2)], 


) 

=S 
1  A 

0<  ir  +  0E; 

(m) 

1>BA  = 

ft  <8^  . 

f>e 

>  -(ir  +  Qg)  ; 

D(m)  - 
wA  - 

CWA  ^VB(pE*  "0E+0, 

0Ei6-2ir“0E  ; 

D(m,= 

WB 

CWA^  [VB(PE' 

2*-0E_e  n 

“0E  —  0—  ■(2»-9g)  ; 

(m) 

°SA  « 

ft1!  tvB(pE. 

V0’  ns)] 

(m) 

(m)  . 

"  "8e1  el0E  ; 

DSB  = 

dsa  (e=-8)  , 

J 

ft- 

C(AT1!>vB.P,  | 

-K8  +9,  2)  + 
E 

vB(pl,~-{i+2)eE+0f2)], 

-  -j  <  0<  w  +  0g  ; 

D{„m)  =  D(m)  (e=-6}t  ~>  Q>  -  (t  +  %)  ; 

Bl  Ai  2-  - 

Diw  ■  ‘ft*  ’+  ve- 2,!-  -f  1 85  ^ 

4*'  ■  ft  (e=-e>-  |>e>  -<»+%>; 
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(17) 

cont. 


=  C(“”1>  [v^p*.,  tt  -0r+O,2)j  ,  -  \  <  G<  tt  +  0e  ; 

AO  ^  t  "" 


and 


“SS’-C*"* 


|>  e>  -(ir  +  eg)  , 


where  the  property  of  symmetry  with  respect  to  9  =  0  is  used  to  obtain 
symmetrical  terms.  The  coupling  coefficients  of  different  orders  are 
obtained  in  Appendix  C.  As  soon  as  coupling  coefficients  are  properly 
evaluated,  the  diffraction  of  any  order  can  then  be  obtained  by  making 
use  of  Eqs.  {16}  and  (17). 


Finally,  the  total  far-field  pattern  of  the  corner  reflector  can  now 
be  obtained  by  superposition  of  all  terms  presented  above.  Taking 
wedge  A  as  a  common  phase  reference,  and  considering  only  the  upper- 
half  region,  0<  0^ir,  the  total  far-field  u(6)  can  be  written  in  a  form 
similar  to  Eq.  (46)  as 


{18} 


u{8) 


v*+  yDr’ 

yAS  + 

_  — 

V  Jm) 

Z  °A 

+ 

f  Jm) 
>  db 

tz  J 

_m=l 

_m=l 

yAB 


0* 

h 

Ot 

V  J™) 

V 

+ 

l  °w 

yAW  +  L 

Z  <Il  1 

jn=2 

i=l 

_in=l 

l  ttu  >h 


m~l 


yBh  yAB, 


where  the  last  images  in  the  upper  wall  arc  included  because  in  general 
they  may  contribute  to  the  upper-half  region.  The  local  phase  factors 
referred  to  A  can  be  written  from  Figs.  1  and  3  as 


yAS  =  ExP*  H2*Pe  cost-Sg+e)]  , 
=  Exp.  I  -j2ir  b  sin  6], 


(19) 


Yaw  =  ^P*  [-j2wpE  cos  (%-0>j  =  YAS> 
y Xi  -  Exp*  L -j2irpj  sin  (i&E+O)], 

and 

yBh  =  Exp.  [  -j2irp.  sin  (hOg-G)]  . 


IV.  THE  APPROXIMATED  SOLUTIONS 


The  pattern  of  the  corner  reflector  in  Fig.  1  can  be  ca'.  -tilated 
by  Eqs.  (18)  and  (19)  a*  accurately  as  desired.  Since  the  contributions 
of  the  higher-order  terms  to  the  pattern  decrease  with  increasing  order, 
computation  can  be  made  by  including  only  those  terms  which  are  signi¬ 
ficant  in  their  defined  regions.  The  following  approximations  are  made 
to  obtain  a  pattern  including  only  significant  higher-order  terms. 


First,  referring  to  Fig.  2a,  the  contributions  from  wedge  S  are 
the  geometrical  optics  term  v  of  Eq.  (1)  due  to  the  primary  source, 
and  the  diffracted  terms  Dg11  due  to  the  induced  sources.  The  lowest- 
order  diffracted  terms  are  due  to  the  second-order-induced  sources. 


Therefore,  the  most  significant  contribution  from  the  diffracted  terms 
is^^jg,  given  in  Eq.  (14).  Computation  of  the  illumination  intensity, 

CgA,  indicates  that  for  typical  dimensions  the  magnitude  of  Dgis 
generally  very  small  compared  to  the  unit  intensity  of  v  .  Neglecting 
the  diffracted  terms,  the  contribution  from  wedge  S  may  be  approximated 
as 


(20) 


* 


v 


0* 

*  y 

L~ 

m=2 


(m) 

% 


«N#  V*  =  1, 


0  <  8<  9]?.* 


where  the  region  of  contribution  is  defined  as  such  because  only  the 
upper-half  region  0<  ir  is  under  consideration. 

Next,  consider  the  con*  ributions  from  the  induced  sources  at  A. 

Hie  third-  and  higher-order  i^.-ms  of  a  typical  reflector  are  also 
generally  very  small  compared  to  unity.  Therefore,  the  diffpacted  fields 
from  A  can  be  approximated  by  considering  only  D^.  and  D in  Eqs.  (2) 
and  (14),  respectively.  The  diffracted  fields  from  B  can  be  approximated 
in  the  same  manner.  We  thus  have 
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0  <  0  C  ir , 


a* 

h-1 

(21) 

> 

n  m)  - 

da  ^ 

D{1) 

S 

•4- 

+  DAB 

V' 

+  > 
L.- 

d(*) 
Ai  * 

m*l 

i=l 

and 

c* 

h-1 

v 

/ 

tw 

(n  ) 

db  % 

D(1) 

VBS 

+  Dfe) 

+  BA 

DU> 
uBi  ’ 

m=l 

i=i 

n 

0<8<  j 


MTr+6E)lei  * 


The  lowest-order  contributions  from  the  induced  sources  at  wedges 

hich  are  also  very  small  compared  to  the  unit 
defined  in  the  same  region  as  v  ,  the  con- 


W  and  S  are  D^)  and 


m: 


D§is 


intensity  of  v  .  Since 

tributton  from  Dg2\o  the  total  pattern  can  be  neglected.  Buts 
defined  outside  the  region  of  v  ;  therefore,  in  order  to  accurately  cal 
culate  the  side  and  back  radiations  tn  the  pattern,  we  must  at  least 
approximate  the  contribution  from  W  as 


(22) 


os 

I 


(m)  ,  . 

Dw  *  °WA  +  ^WB 


,(*) 


$£<  6<  ff 


m=2 


As  for  contributions  from  the  image  terms,  the  magnitude,  in 
general,  decreases  with  increasing  value  of  i.  In  other  words,  in  the 
region  of  contribution,  0<  8<  ttlZ  -  9g,  shown  ir.  Fig.  3,  the  images 
contribute  to  the  total  pattern  more  significantly  in  the  upper  end  than 
in  the  lower  end  of  the  region.  Therefore,  in  theory,  approximation  on 
the  images  should  be  made  individually.  But,  we  shall  approximate 
images  of  the  same  order  as  a  group  so  that  all  the  images  of  the  same 
order  can  be  completely  included.  Assuming  that  the  second-  and  higher- 
order  images  have  negligible  contribution  to  the  total  pattern,  the  image 
term3  are  approximated  as 


m 


|  -  -  iDj. 
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An  approximate  pattern  of  the  corner  reflector  can  now  be  ob¬ 
tained  from  Eq.  (18)  by  using  Eqs.  (19)  to  (23),  as  follows: 


u{9)  £  l 


v*iyAS  +  !Das  +  °AB  +  J  1 


+  [dbs  +  dLa+  f  db!  I  >'AB 


+  [D^  +  1  y 

1  WA  WB ' yAW 


41 

+  1  UL)ili'Al+[(IU1)hlyBh>rJ 


Since  this  is  an  approximated  pattern,  discontinuities  are  expected  to 
be  increasingly  noticeable  with  decreasing  size  of  the  corner  reflector. 
Let  ua  examine,  term  by  term,  the  continuity  of  Eq.  (24)  in  the  upper- 
half  region,  0  £  0£  tr.  At  6s  the  discontinuity  of  v  in  Eq.  (20)  is 
eliminated  by  «  At  8  =  ir/2,  the  discontinuity  of  in  Eq.  (21)  is 

eliminated  by  ;  but  there  are  no  higher-order  terms  included  in 

Eq.  (24)  to  compensate  for  the  discontinuities  of  D^\  and  D1;,  with 

BA  Bi 

*  =  1,2,  —  (h-1).  Similarly,  at  9  =  -(?r  +  9jr),  the  discontinuity  of 

is  eliminated  by  but  those  of  and  D^,  with  Is  1,2, - (h-1), 

are  left  uncompensated.  The  discontinuities  at  9  =  ir/2  -  i6g  of 

are  all  eliminated  by  and  ,  with  i*  1,2,  --  -(h-1).  The  dis- 
{li  At  Bi  * 

continuity  of  (1^  )h  takes  place  in  the  defined  regions  of  v  and  is  usually 

unnoticeably  small.  The  discontinuities  of  and  D^.g  at  9  =  8g  are 

also  unnoticeably  small  for  a  typical  corner  reflector. 

Within  the  accuracy  of  the  approximations  made  to  obtain  Eq.  (24), 
i-  is  desirable  to  have  the  pattern  continuous  in  the  entire  region.  To 
accomplish  this,  the  second-order  discontinuities  at  9  =  x/2  euid  -(x  +  %) 
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mentioned  above  need  to  be  eliminated.  The  coupling  coefficients  from 
E  to  A  resulting  from  the  second-order-induced  sources  at  B  can  be 
obtained  from  Eq.  (14)  as 


,25)  CW.D»(6  =  |). 

=  CBA  lvB<b’  2)  +  vB(b’  *  '  2b£‘ 
h-1 

+  )  ^Bi  ^  vB^i’  i0E'  +  vB^i*  v~  (i+2)0£»  2)}  , 
i=l 

h-1 

—  c  +  'S'  c 

-  CBA  CABA  +  /  CBi  CABi 

i-1 


-  c(1)  C  + 
"  CAB  CABA 


h-1 

z 

i=l 


*0 


C' 

At 


'ABi 


where  Eqs.  (9)  and  (10)  are  used  and  property  of  symmetry  is  employed 
for  the  first-order  coupling  coefficients.  By  symmetry,  the  second- 
order  coupling  coefficient  from  A  to  B  can  also  be  obtained  from  Eqs . 
(14),  (7),  and  { 1 0)  as 


(9  = 


(*> 

AB* 


Because  of  the  interaction  between  A  and  B  by  C^g,  a  third-order  dif¬ 
fraction  can  be  written  for  A  and  B  with  the  same  defined  region  in  Eq- 
(14).  If  the  process  of  interaction  between  A  and  B  continues  to  infinite 
order,  a  coupling  coefficient  can  be  obtained  in  closed  form  as 


c<»>  +  ¥ 
AE+  L 

i=l 


(26) 


'AB 


r(i) 

°Ai 


'ABi 


1  -  C 


ABA 


*9 


where  Ca3A  and  CABi  can  obtained  from  Eq.  (25).  The  continuity 
of  the  total  pattern  at  8*  ir/2  is  now  ensured  by  using  C^g  instead  of 


for  D<0  ar.d  D{*\  ™  Eq.  (24) 


AB 


BA 


as 


(27) 

and 


DAB  *  CAB[vB(b,  y  +  6,2)  4  vB(b,-^  -  28E+6,  2>],  °£-9£». 


DBA  *  CAB  1  -6’  21  +  vBlb>  ¥  -26E-6.  2!).  | 


0  <  9sL 
-(ir+6g}<  G<  n , 


where  £q.  (?)  is  used  and  the  regions  are  restricted  to  the  upper-half 
region. 

As  a  consequence  of  the  modified  equations  in  Eq.  (27),  the  un¬ 
compensated  term3  at  0  =  -(lr  +  Op.)  are  now  D^.  and  with  i  1  *  2 

15A  Bi  # 

-  —  (h—  1 ) -  The  coupling  coefficient  from  B  to  W  resulting  from  these 
terms  can  be  written  from  Eqs.  (27)  and  (10)  as 


cTO  =  DBA(e=-("+9E»+  >  ngj  <»-<••+ He)) 


h-1 


CWA=DAB,e=  ”  +  eK>+  I  DA!f0=,r+eE> 


i=l 


*  CAB  ¥  +  0E'  2>  +  ~T  -  9E'  2>J 


'AB 


h-1 

+  J  Ca!  [vB(Pj,  ¥  +  (I+1>SE-  2)  +  vB(Pf.  ¥~  ii+1)0E-  2))  > 
i=l 


which  is  a  new  coupling  coefficient  to  ensure  the  continuity  of  the  total 
pattern  at  9  =  -(-ir-fOg)  or  ir-9g.  Adding  this  new  coupling  coefficient  to 

^"WA  *n  ^**1*  (7),  the  diffracted  fields  at  W  are  modified  as 
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M  * 


(28)  %a  =  ^WA^v3^E'  “eE  +  8*  nW^ 


D, 


where 


WB  =  CWA^Vb^’  2ir -0p,-8,  nw)j 


CWA  =  C^A  +  C^A 


Using  Eqs .  (27)  and  (28),  an  approximated  "continuous"  pattern 
can  be  finally  written  from  Eq.  (24)  as 


h-1 

(29)  u(0)  =  [v*)yAS  +  (D^  +  Dab  +  I  DAi  ‘ 


+  [Dbs  +  dba  + 


h-1 


Y  DM)y 

/  ^BiiyAB 


i=l 


+  [Daw  +  ®bw1  yaw 


h 

*  Y  ^*L^yAi  +  [^U  MyBbyAB’ 

i=l 

For  typical  dimensions  of  a  corner  reflector,  Eq.  (291  in  general  gives 
excellent  prediction  of  the  radiation  pattern.  Examples  are  given  in  the 
following  section. 


V.  COMPUTED  PATTERNS  COMPARED  WITH 

MEASURED  PATTERNS  OF  HORN  ANTENNAS 

To  illustrate  the  validity  of  the  corner  reflector  as  a  model  of  the 
pyramidal  horn  antenna  fed  by  a  waveguide  supporting  TEio  mode,  Eq. 
(29)  is  computed  and  compared  with  measured  patterns.  Figure  4  shows 
the  experimental  set-up  of  a  horn  antenna  in  which  the  idealized  model 
is  the  corner  reflector  ASB  used  to  derive  Eq.  (29).  The  associated 
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Fig.  4,  The  E-plane  of  a  horn  antenna. 


waveguide  and  the  calibrated  attenuator  and  diode  detector  are  not  con¬ 
sidered  in  pattern  prediction. 

First,  consider  a  horn  antenna  of  pg.  =  41 .3  cn*  an°  =  35“  fed 
by  a  waveguide  propagating  the  TEio  mode  at  9.8  gc.  The  horn  length, 

Pj.,  in  terms  of  wavelength  is  equal  to  13.5X  .  The  measured  far-field 
pattern  is  shown  in  Fig.  5.  The  pattern  computed  by  Eq.(29)  is  shown 
displaced  5  db  below  the  measured  pattern.  Comparison  of  two  patterns 
shows  excellent  agreement  in  the  overall  lobe  structure.  The  small 
deviation  of  relative  field  intensity  in  the  region  50“  <  6<  80“  is  primarily 
due  to  the  approximation  assumed  in  Eq.  (23)  that  the  second-  and  higher- 
order  images  are  negligible.  The  presence  of  the  waveguide  and  the 
associated  attenuator  and  detector  shown  in  Fig.  4  is  responsible  for  the 
interference  in  the  region  80“  <  8<  18CP  of  the  measured  pattern.  In  the 
process  of  computation,  special  care  is  required  for  the  properly  de¬ 
fined  regions  of  each  term  in  Eq.  (29).  Attention  is  also  needed  for 
the  boundaries  S=  dj-,  90*,  (90*-i9j.),  and  (180“-8g.)  at  which  discontinu¬ 
ities  are  eliminated  by  higher-order  terms.  For  convenience,  the  pat¬ 
tern  computed  by  only  first-order  diffraction  treated  in  Ref.  5  is  plotted 
displaced  5  db  above  the  measured  pattern.  As  mentioned  earlier,  the 
discontinuities  at  9  =  90“  and  (180*-8£)are  expected  in  the  first-order 
pattern. 
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A  comparison  of  the  patterns 


When  the  frequency  is  increased  corresponding  to  p  g.  =  24. 8 X 
for  the  same  horn*  the  three  patterns  are  as  shown  in  Fig.  6.  The 
same  conclusions  drawn  for  Fig.  5  remain  true,  except  that  the  inter¬ 
ference  from  the  waveguide  and  the  associated  structure  tecomes 
larger  because  the  physical  size  is  larger  in  terms  of  wavelength. 

In  Fig.  7 ,  three  patterns  are  shown  for  a  small  horn  antenna  of 
»E  =  5  -  6 IX  and  28g.  =  21.2*.  Although  the  overall  lobe  structure  is 
still  In  good  agreement  with  the  measured  pattern  ,  a  larger  deviation 
in  intensity  level  is  observed  around  8  =  80#  of  the  pattern  by  Eq.  (29). 
This  disagreement  results  because  the  second-order  image  terms 
neglected  in  Eq.  (23)  are  not  negligibly  small  for  small  horns.  There¬ 
fore,  better  patterns  can  be  obtained  for  small  horns  by  including  the 
second-order  image  terms  and  their  subsequent  effects  on  the  total 
pattern  . 

The  three  examples  presented  above  have  demonstrated  the 
accuracy  of  Eq.  (29)  for  pattern  computations  of  typical  horn  antennas. 
The  accuracy  of  she  experimental  measurements  is  assumed  to  have  1  db 
fluctuation  when  the  intensity  is  around  40  db  below  the  reference  inten¬ 
sity.  In  view  of  this,  Eq.  (29)  is  sufficient  for  horn  antennas  of  typical 
dimensions .  When  pj-  and  6g-.  become  smaller,  it  is  easily  observed  from 
Fig.  7  that  the  second-order  image  terms  in  Eq.  (15)  should  be  in¬ 
cluded  to  ensure  good  prediction  around  the  region  8i  90“  -  8g. 


VI.  RELATIVE  BACK  LEVELS 

The  radiation  patterns,  either  measured  or  computed,  always 
have  a  back-lobe  maximum  at  9  =  180",  even  though  this  maximum  value 
is  not  necessarily  the  largest  maximum  in  the  region  ISO*  -  8g  <  8<  180*. 
In  Fig.  6,  for  instance,  the  largest  maximum,  value  in  the  region  is  at 
176*  of  the  measured  pattern  and  at  8  =  172.5*  of  the  computed  patterns, 
using  Eq.  (29).  This  discrepancy  results  from  the  interference  from  the 
experimental  set  up  which  is  not  considered  in  the  idealized  reflector 
model.  The  difference  between  the  value  at  9  =  180*  and  that  of  the 
largest  maximum  is  generally  small.  Therefore,  the  radiation  intensity 
at  3  =  180*  can  be  taken  as  a  representative  value  for  back-lobe  region. 

It  is  also  observed  that  the  radiation  intensity  in  the  region  90*  <  8<  180* 
-8g  is,  in  general,  smaller  than  the  mentioned  representative  value. 

In  EMC  problems,  it  may  be  desirable  t,o  keep  the  pattern  level  in 
the  region  90*  <  0<  180*  as  low  as  possible.  To  predict  a  representative 
interference-intensity  in  this  region,  one  can  simply  take  the  first-order 
diffraction  terms  in  Eq.  (2)  for  approximate  computation.  Neglecting  all 
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Comparison  of  patterns 


Comparison  of  patterns 


higher-order  terms*  the  radiation  intensity  at  0  =  ISO*  is  given  approxi¬ 
mately  by 


u(")  =  (0=ir)  +  (0=  -ir) 


=  2  vg(pE»  2n  -  ©E»  2)  . 


The  radiation  intensity  at  0  =  0*  is  also  approximated  as 

u(0)  =  Xas(0=°)  +  D^(e=°)  +  £>gs  (9  =  0) 

=  Exp[  -j2ir  p£  cos  Og]  +  2  vB(pE,  n  -  0E,  2). 


where  Eqs.  (19)  and  (2)  are  used.  The  back-to-front  ratio  can  then  be 
obtained  by  plotting 


20  log 


)u(*) 1 

|u{0)J 


which  has  only  two  parameters,  0£  and  pE.  For  a  horn  antenna  of  fixed 
angle  0g,  the  back-to-front  ratio  can  be  plotted  as  a  function  Figure 

8  shows  seven  curves  for  different  horn  dimensions  to  demonstrate  that 
Fig.  8  can  be  used  for  approximate  prediction  of  typical  horn  antennas, 
Fig.  9  gives  a  set  of  measured  data  for  20E=  35* ,  The  data  are  taken 
from  nine  complete  patterns  of  tne  antenna.  The  accuracy  of  prediction 
depends  on  both  the  approximation  in  computation  and  the  equipment  used 
in  measurement.  In  general,,  it  is  expected  that  the  larger  the  sises  of 
antenna  the  better  tne  prediction. 


It  is  interesting  to  compare  the  first  minimum  points  in  Fig.  8  to 
the  defined  optimumf  6]  horn  lengths.  The- values  of  optimum  horn  lengths 
are  tabulated  in  Fig.  8  and  are  found  to  be  one  half  of  that  of  pg,  which 
give  the  first  minimum  back-to-front  ratios . 

The  sets  of  pE  and  0£  giving  rise  to  the  minima  in  the  curves  do 
not  necessarily  imply  that  u(0)  is  maximum  and  that  u{ir)  is  minimum. 

The  reason  for  this  is  the  main  lobe  of  pattern  begins  to  bifurcate  at  the 
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mum  ) 


Fig,  9.  Back-to-front  ratio  of  2  0£  =  35*. 


points  where  minimum  back-to-front  ratios  occur.  For  a  horn  antenna 
with  90  =  25*,  and  p£  =  6.  5k,  17k,  27k,  38k,  and  48k,  where  minima 
take  place  in  Fig.  8,  the  main  beam  of  the  pattern  is  split  into  two, 
four,  six,  eight,  and  ten  looes,  respectively.  As  for  EMC  problems. 
Fig.  6  can  be  used  to  evaluate  approximately  the  representative  back- 
radiation  intensity  u{ir)  relative  to  front- radiation  intensity  u{Q),  If 
the  horn  antenna  is  large  enough,  it  ia  generally  safe  to  expect  that 
the  radiation  intensity  on  the  average  in  the  region  90*  <  9<  180*  - 
is  about  5  db  lower  than  the  representative  value  u{w)  as  shown  in  Figs. 
8  and  9. 


VH.  CONCLUSIONS 

The  E-plane  patterns,  including  far-sidelobes  and  backlobes,  of 
horn  antennas  has  been  formulated  without  employing  aperture  method# 
and  equivalence  principles.  Considering  the  various  assumptions  and 
the  mathematical  difficulties  inherent  in  aperture  methods  it  is  shown 
here  that  diffraction  theory  is  more  accurate  and  practical  in  analysing 
radiation  characteristics  of  typical  horn  antennas. 

The  pattern  of  the  proposed  reflector  model  may*  in  theory,  be 
computed  by  Eq,  (18)  as  accurately  as  desired.  For  reflectors  of 
typical  dimensions,  the  approximated  pattern  in  Eq,  (29)  has  been  shown 
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in  i'igs.  5,6,  and  7  to  be  in  excellent  agreement  with  the  measured 
patterns  of  norn  antennas  for  which  the  reflector  model  is  intended. 
Comparisons  with  the  first-order  patterns  indicate  that  the  improve¬ 
ments  of  Eq.  (29)  are  mairdy  in  the  far-sidelobes  and  backlobes.  As 
a  consequence  of  the  approximations  made  to  obtain  Eq.  (29),  the  pat¬ 
tern  level  tends  to  deviate  more  and  more  in  the  region  around  9=  90* 
-0E>  when  horn  dimensions  become  smaller.  Figure  7  shows  that  the 
second-order  image  terms  in  Eq,  (15)  can  no  longer  be  assumed  negli¬ 
gible  for  small  horns.  In  conclusion,  Eq.  (29)  is  generally  sufficient 
to  predict  patterns  of  typical  horn  antennas.  If  higher  accuracy  is 
desired,  the  higher-order  terms  neglected  in  Eq.  (29)  may  be  included 
for  computation. 

The  curves  of  backlobe  levels  in  Figs.  8  and  9  are  approximated 
by  the  first-order  diffraction  terms.  More  accurate  prediction  of  back- 
to-front  ratios  should  include  all  higher-order  terms  which  have  signifi¬ 
cant  contribution  in  the  back-lobe  region. 
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APPENDIX  A 

OF  DIFFRACTION  THEORY 


The  two-dimensional  problem  of  the  electromagnetic  field  in 
the  neighborhood  of  a  conducting  wedge  illuminated  by  a  uniform  plane 
wave  was  first  solved  by  Sommerfeld(7 ]  •  The  solution  for  a  conducting 
half-plane  (zero  wedge  angle)  was  formulated  in  terms  of  the  Fresnel 
integral.  Subsequently,  Pauli[8]  formulated  the  solution  for  wedges 
of  arbitrary  angles  in  an  asymptotic  series  in  which  the  dominant  term 
is  the  Fresnel  integral.  The  higher-order  terms  in  Pauli’s  solution 
become  identically  zero  for  zero  wedge  angle.  Therefore,  Pauli's 
solution  is  used  here  for  the  general  case. 

Figure  10a  show’s  the  geometry  of  the  wedge  used  by  Pauli  to 
formulate  the  solution  of  field  intensity  at  P{r,  40  caused  by  plane-wave 
illumination.  By  reciprocity,  if  the  same  wedge  of  a  perfectly  conduct¬ 
ing  surface  is  illuminated  by  a  uniform  cylindrical  wave  from  S  shown 
in  Fig.  iOb,  the  far-field  intensity  can  be  written  from  Pauli,s  solution 
as 


(30) 


v  =  v(p,  <P  .  n)  +  v(p,  4>‘.  n), 


=  v+  i  v 


and 


TT 


where  the  terms  represent  the  incident  and  the  reflected  fields, 
respectively.  The  sum  (v*  +  v“)  applies  when  S  is  a  magnetic  line  source, 
and  the  difference  (v+  -  v“)  applies  if  S  is  an  electric  line  source.  The 
incident  and  the  reflected  fields  are  composed  of  geometrical  optics 
terms  and  diffracted  terms,  as. 


v±  =  (v*)*  +  v*  . 
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(a)  Illuminated  by  a  uniform  plane  wave 


(b)  Illuminated  by  uniform  cylindrical  wave 
Fig.  10.  Geometry  of  a  wedge. 


S  >’  *  1 1, 


The  geometrical  optics  terms  are  given  by 


(31) 


(v*)*  =  v*{p,4>^,  n) 


Exp[  jkp  cos(<j>*  +  2x nN)]  »  -ir  <  +  2wnN  <  *•, 

0  otherwise 


N  =  0,  *1,  ±2 


where  the  periodicity  of  the  functions  is  seen  to  be  2nw.  The  diffracted 
terms  are  given  by 


(32) 


VB*  =  vb(P*  n) 


2  j  cos-y-  j  Exp  (jkp  cos  $*) 


2 


dr 


cos 


cos- 


4> 


+ 


(Higher-order  terms 
negligible  for  large  kp 


and 


a  =  1  +  cos  <f>- 


As  (akp  )  approaches  infinity,  i,  e. ,  as  the  line  source  S  recedes  to  the 
far-zone,  the  solutions  in  Eq.  (32)  can  be  written 


(33) 


vB(p,$  ,u)  = 


Exp[-j(  f  +kp)} 
,j  2irkp 


sin  ~ 
n 


i  *  A-  . 

n(cos  —  -  cos  a —  ) 
n  n 


For  N=0  in  Eq.  (31),  the  geometrical  optics  terms  have  discontin¬ 
uities  at  the  shadow  boundary  (<j>+  =  *}»  and  the  reflection  boundary  (<j>”  **). 
At  these  boundaries,  the  diffracted  terms  are  given  as 
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(34) 


vB(*  0*=tt,  n)  = 


i 

-  J  Exp.  {-jkp}  + - 

\  i 

+  7  Exp.  (-jkp)  + - 


<j>*  ~  lim{TT  -t) 

t  —0 

6^  =  lim  (w  +c ) 
c  -0 


where  the  higher-order  terms  are  negligible  for  large  (kp)  are  not 
presented.  The  solutions  are  essential  in  ensuring  continuity  at  the 
boundaries.  Since  we  intentionally  make  the  geometrical  optics  terms 
in  Eq.  (31)  defined  at  the  boundaries,  =  +ir,  the  first  equation  of 
Eq.  (34)  should  be  used  to  obtain  the  field  intensities  as.  the  boundaries. 
The  total  far-field  of  the  wedge  illuminated  by  a  uniform  cylindrical 
wave  shown  in  Fig.  10b  can  now  be  obtained  by  using  Eqs.  (31),  (32), 
and  (34)  as 


(35)  u(p,  <p,  n)  a  [(v*)+  +  vB+]  i?4(v*)-  A  vg  J  , 


where  the  reflection  terms  disappear  if  the  wedge  is  illuminated  by  the 
source  with  ip  =  0.  If  the  source  is  an  electric  line  element,  the  value 
of  u(p,  4>,  n)  is  identically  zero  for  =  0. 

The  solutions  in  Eq.  (35)  are  valid  only  in  evaluating  the  far-field 
intensity  which  is  the  main  concern  of  present  problem.  The  diffracted 
near-field  intensities  of  a  conducting  half-plane,  i.  e. ,  n  =  2,  have  been 
solved  by  Nomura[9l  and  used  by  Ohba  for  dipole  source  illumination. 
Since  this  solution  is  in  a  general  form,  it  can  easily  serve  to  illustrate 
principles  of  reciprocity.  To  generalize  the  solution  further,  the  dif¬ 
fracted  near-fields  of  a  wedge  illuminated  by  a  line  source  have  been 
written  by  Dybdal(l0|  in  a  form  as  follows: 


(36)  vlr,,. <.*,„)=  SfckggEMl  [VB(h>  ♦*.„)] 

r  +  p 


JLEL 

r+p  ’ 


where  the  vg  terms  art  given  in  Eq.  (32).  If  the  field  point  is  in  the 
far- zone,  v(r,  p,  4>*,  n)  is  reduced  to 


34 


v(r,  p,  n)  *  I^P.  [  vB(p»  <+>“,  n)j  , 

\ !r 

i 

which,  with  the  common  factor  r" *  Exp(-jkr]  removed,  is  identical  to 
Eq.  (32),  If  the  source  point  recedes  to  the  far-zone,  the  solutions  in 
Eq.  (36)  are  reduced  to 

v(r,p.  <►*„>*  &2L&£l  [vB(*,«>*.n)]  , 

|P 

which,  with  the  common  factor  p~*  Exp[-jkp]  removed,  are  identical 
to  the  near-field  solutions  of  a  wedge  illuminated  by  a  uniform  plane 
wave. 
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APPENDIX  B 

FAR-FIELD  WAVES  FROM  THE  IMAGES 
OF  LINE  SOURCES 


Consider  line  sources  parallel  to  perfectly  conducting  plane 
walls.  The  problem  in  this  Appendix  is  to  obtain  the  far-field  intensity 
of  the  image  waves  from  ‘he  walls.  The  problem  is  a  two-dimensional 
one  and  the  line  sources  can  be  either  magnetic  or  electric  line  sources. 

First,  a  directional  line  source  is  shown  at  point  S  in  Fig.  H, 

If  a  semi-infinite  conducting  wall  is  placed  with  an  angle  8Q  with  respect 
to  the  horizontal  reference  axis,  the  image  of  the  line  source  is  formed 
at  point  1.  The  line  source  radiates  a  directional  cylindrical  wave  to  the 


HALF- PLANE  WALL 

Fig.  11.  Geometry  of  a  line  source  and  its  image. 


far-field  and  it  is  designated  as  v(8).  In  the  presence  of  the  wall,  v{8) 
is  reflected  by  the  wall  in  the  region  -{ir  +  8Q)  <  8<  -  it/2.  The  reflected 
wave  is  called  the  image  wave  from  point  1  and  can  be  described  by 

(37)  L  (®)  “  v“(6)  ~  v(-26o-0)  ,  §  -  28Q  <  8<  -(*  -  0O)  , 

where  the  superscript  minus  sign  implies  the  reflection  of  v{9).  The 
total  far-field  pattern  can  then  be  obtained  by  superposition  of  the  two 
waves  as 
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(38)  u(9)  =  v(*)  +  v(-20o  -  0)  y8i 

and 

ysl  =  Expl-jkp  sin  (0+  0O)J  , 


where  ysi  is  the  local  phase -factor  of  the  image  referred  to  point  S 
and  p  is  the  distance  between  the  image  and  the  source  as  shown  in 
Fi^.  11.  This  result  is  obtained  by  removing  the  common  factor 
R-a  Exp  [  -jkR] ,  with  R  referred  to  S.  The  plus  sign  applies  for  a 
magnetic  line  source,  while  the  minus  sign  should  be  used  for  an 
electric  line  source. 

Next,  a  corner  reflector  of  20Q  is  shown  in  Fig.  12a.  Let  there 
be  only  one  line  source  v(8)  at  A.  In  the  region  -(*  -  0O)  <  6<  -w72>  the 
rays  from  A  are  reflected  by  the  lower  wall.  Only  rays  in  the  ray  -1 
zone  (from  image  -1)  are  directly  reflected  to  the  far-field.  The  rays 
in  ray  -2  zone  are  reflected  twice  while  those  in  ray  -3  zone  are  re¬ 
flected  three  times  inside  the  reflector.  Consequently,  three  image 
waves  are  formed  in  three  distinct  regions.  The  first  image  wave  can 
be  written  similar  to  Eq.  (37)  from  Fig.  12a  as 


(39)  vi‘(0)  =  v(-2eo-0)  ,  §-2eo<0<f-eo, 

which  is  obtained  by  replacing  the  0  of  v(B)  by  (-20o  -  0).  The  second 
image  is  formed  in  the  xipper  wall.  This  image  wave  can  be  obtained 
by  replacing  the  B  of  v*"(S)  by  (+20o  -  0)  as 


(40)  vz-(0)  =  v(-40rJ  +  6),  -  (~  -  30o)  9>  -  -  20o)  . 

Again  replacing  the  9  of  the  above  expression  by  <-20o  -  0),  the  third 
image  wave  can  be  obtained  as 

v,‘(0)  =  v(-60o  -  0),  §  -  40o  <  6<  *  -  780  . 


This  process  can  be  used  for  any  number  of  images.  The  number  of 
Images  is  equal  to  the  number  of  ray  zones  determined  by  the  highest 
integer  h  such  that  h  <  *72 0O.  Construction  of  the  images  in  Fig.  12a  can 
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be  extended  to  conclude  that  all  the  odd-numbered  images  are  formed 
in  the  lower  wall,  while  the  even-numbered  images  are  in  the  upper 
wall.  In  general,  the  image  waves  can  be  written 

v{-2i0o  _  0)  --  -(i+1  )0Q  <  G<  —  -  i0Q  for  i  odd 

(41)  vi-(0)  =  L  L 

v(-2i0o  +  0)  _[j-(i+l)0o]  >  8>  -  [  |--i0o]  for  i  even, 

i  =  1,2,3 - (h-1)  . 

If  the  ratio  of  x  to  20o  is  not  exactly  an  integer,  the  valid  region  of  the 
last  image  should  be  modified  as 

v(-2hGo-0)  7  -(h+l)0o<  G<  ir  -{2h+l)0o  h  odd 

(42)  vh'(?)  = 

v(-2h0o+0)  -[  1  -(h+l)0o]  >  0>  - [tt  -(2h+l)0o]  h  even  . 

_  2 

The  value  of  20o  in  Fig.  12  is  larger  than  v/4  but  smaller  than  vlZ. 
Therefore  h  is  equal  to  3  and  the  last  image  is  defined  by  £q.  (4:2) 
for  h  odd. 

Figure  12b  shows  two  symmetrical  line  sources  at  A  and  B  from 
which  far-field  waves  v(0)  and  v(-0),  respectively,  are  radiated.  The 
images  of  the  source  at  A  are  identical  to  those  in  Fig.  12.  Because 
of  the  symmetrical  properties  of  the  assumed  line  sources  and  the  geo¬ 
metry  of  the  reflector,  the  images  resulting  from  the  source  at  B  are 
symmetrical  to  those  of  A.  Consequently,  equal  numbers  of  symmetri¬ 
cal  images  are  formed  in  both  walls  of  the  reflector.  Making  use  of  this 
symmetry,  the  image  waves  excited  by  v(-0)  from  B  can  be  obtained 
similar  to  Eq.  (41)  as 

[|-(i+i)0j  >  e>-[I-ie0]  iodd 

2  -(i+l)e0  <  0<  7  -i0o  1  even. 


In  case  the  ratio  tt  to  20o  is  not  an  integer,  the  last  image  wave  has 
its  defined  region  as  in  Eq.  (42)  by  interchanging  h  odd  and  h  even. 
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Comparison  of  Eos,  {41}  to  {43}  indicates  that  all  image  waves 
from  the  lower  wall  are  combinations  of  v{(d)  with  i  odd  and  v{ (-9) 
with  i  even.  The  other  set  of  combinations  gives  the  image  waves 
from  the  upper  wall.  To  summarize,  the  image  waves  from  both  walls 
of  the  reflecto:  ;ir  "'esignated  by  <IjJi  and  (lyU  as  follows: 

(44)  (IL)i  =  v“(9)  =  v(-2ieo-6)  ,  |~(i+l)0o<  0l|-i9O'  lower  wall; 

(Iu)i  =  v"(_8)  =  v{-2i0o+0)  .  -[  |  ~{t41)0o]  >  «Ppe*  wall 

i  =  1,2,3, - h  . 

For  the  ratio  of  n  to  20o  not  an  integer,  modifications  for  i  =  h  are 

(45)  (lL)h  =  v(-2h0o-6),  l-{h+l)0o<  0<  ir>(2h+l)0o; 

and 

(Iu)h  =  v(-2h0o+0)t  -[I.{h+l)8e]  >  0>  -[ir-(2h+l)0oj  . 

w 

In  Fig.  I2b,  three  images  in  the  lower  wall  are  shown.  Their  far-field 
intensities  can  be  obtained  by  setting  h  =  3  in  Eqs.  (44)  and  (45). 

Because  of  symmetry,  only  the  upper-half  region  0<  0<  needs 
to  be  considered  for  tne  radiation  pattern  of  the  reflector  antenna.  In 
this  region,  the  contributions  are  from  the  sources  at  A  and  B  and  the 
images  in  the  lower  wall.  Contributions  from  the  last  image  in  the 
upper  wall  are  possible,  if  the  ratio  of  w  to  20o  i3  not  an  integer.  Super¬ 
position  of  all  these  contributions  gives  the  total  far-field  u{8)  as 

h 

(46)  u{0)  =  v(0;  +  [v{-«}]y ^  (il)1  [v(-2i0o-0)]yAi 

t=l 

+  (+l)hiv(-2h9o+0)JyAB  yBfc  , 


W 


where  (ii)1  preceding  the  image  terms  are  used  for  magnetic  line  sources 
and  (-1)  are  used  for  electric  line  sources.  The  phase  factors  are 
introduced  by  taking  A  as  phase- reference: 

(47)  YAB  =  ^P  t-jkPo  sin  03* 

yAi  =  ***?  ^  ”ikPi  sin  <i0o+e)3  » 
yBh  =  [’jkPh  sin(h0o-0)J  » 

PQ  =  AB,  and 

Pi  =  Pi_l  cos  0O  +  p0  cos  (i0o). 

The  above  expressions  can  be  obtained  by  considering  the  geometry  of 
Fig.  12,  It  is  noted  that  each  term  in  Eqs.  (4S)  and  (47)  is  set  zero 
outside  its  defined  region. 


APPENDIX  C 

COUPLING  COEFFICIENTS  OF  THE  FOUR  WEDGES 
OF  THE  ANTENNA  MODEL 


The  primary  source  at  S  in  Fig.  1  illuminate*  wedge*  A  and  B. 

The  resulting  diffractions  are  denoted  as  first-order.  When  the  first- 
order-induced  line  source  at  A  illuminates  B,  W,  and  S,  the  diffractions 
are  designated  as  second-order.  Likewise,  A,  W,  and  S  are  illuminated 
by  the  first-order-induced  source  at  B.  Therefore,  the  lowest  order  of 
diffractions  for  W  and  S  in  our  present  problem  is  second-order. 


Physical  significance  of  the  coupling  coefficients  in  Eq.  {17)  can 
be  more  clearly  stated  by  taking  Cy£g“  as  example.  Wedge  B  is 
illuminated  by  the  {m-ZJ^-induced  sources  at  A,  W,  S,  and  (fa-1)  images. 
There  are  then  (h+2)  terms  of  (m-l)^-order  diffraction  waves  induced 
at  wedge  B.  Combination  of  these  (h+2)  terms  forms  a  directional 
cylindrical  wave  radiated  from  B.  In  the  direction  of  wedge  A,  8  m 
+  v/Z,  the  illumination  intensity  is  therefore  a  combination  of  (h+2) 
terms  in  this  direction.  Writing  out  these  terms  yields 


,{m-l) 
AB  * 


(m-2)  +C(m*‘2>C  +C{m"2)C 

BA  '"'ABA  *  ^BW  UABW*^BS  ^ABS 


h-1 

+  ^  CBi °ABi * 
i*l 


where  the  double  subscripts  denote  (m-2)**border  illuminations  from 
various  sources  to  B.  The  triple  subscript  terms  denote  coupling  co¬ 
efficient*  from  various  sources  through  B  to  A.  Similarly,  for  other 
wedge#  and  the  images,  symmetry  property  can  be  used  to  obtain 


(48) 


Cba  l) 


(m-2)  c(R1"2)C  +r(m"2)C 

-  CAB  °BAB  +  °AW  ‘-SAW  +  CAS  CBAS 


+ 


(m-2) 

CAi 


CBAi 


r(m-l) 

CAB 


9 


(48) 

(cont) 


c<m-l)  (m-2)c 
^  Ai  =  Car  '■'Aii  +  ^ 


Cab  CA1A  +  CAW  CAiW  +  CAS  CAiS 
h-i 


Z  Ca^  C 


(m-2) 


Ak  '"Aik  =  BBi 


(m-i) 


k=l 


(rn-I)  (m-2)  {rn- 2}  (m-2) 

AW  =  CWA  AWA  WB  °AWB  ~  BW  ’ 

(m-1)  (m-2)  (m-2)  _  (m-2) 

'-'AS  =  CSA  CASA  ^SB  CASB  “  CBS  * 


'A3 

CWAB  +  CAW 

CWAW 

h-1 

I 

c{m‘2)c 

CAi  CWAi  ‘ 

(m-1) 

CWB 

i=l 

,Jni-2) 


(m-1)  (m-2) 

"SA  3  CAB 


Cca  =  C’.„  '  CSAB  +  CAW  C-  +  CAS  C 


(n»-2)  +  (m-2) 

SAW  AS  ~SAS 


1* 


+  \  (m-2) 

i»l 

m  =  3, 44  5  - » 


CAi  CSAi  ~  CSB 


(m-1) 


Equation  (48)  applies  only  for  m>  3.  For  m  =  1,  the  coupling  coefficient* 
are  of  zero-order  and  are  not  defined  for  our  problem  except  for  direct 
illumination  from  the  primary  source  at  S,  For  m  =  2,  the  coupling  coef¬ 
ficients  are  of  first-order  and  are  given  in  Formulation  of  Solution  as 
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,'t  . . . 


«4»  CBA  =  CZ  -  '»((*  f  ’  9o)  • 


and 


CWA  =  =iJk  *  vb'pe-  2’!  • 
C£a  ~  =  VB*PE‘ 


Ca!  s  CBi  =  vb(pE*  I  -(i+1Jeo)‘  i  *  1.2.--  (h-1)  . 

Thfi  C  and  axe  equal  to  *rro,  for  m  *  2, 

because  there  are  no  first-order  induced  sources  at  W  and  S.  For  m  *  3, 
Eq.  {48}  can  be  written  as 


{50} 


h-1 

a-(s)  _  rW  r  i  Y  r^)  r 

CBA  CAB  °BAB  L  CAt  °- 


1*1 

h-i 


* 

BAi  CAB  * 


-  c ^  C  +  V  c  * 

Ai  ~  CAB  AiB  L  Aik  CAik  CBi  ’ 

k*l 


and 


C{i)  -r{l)  C  +C{1>  C  =C{2) 

AW  “  CWA  AWA  WB  UAWB  BW  ’ 

c<*)  _  _m  r  .  _{!}  _  _  r{*) 

CAS  ~  CSA  CASA  +  CSB  CASB  ’  Ss  * 

h-1 

CW  -  rW  c  +  Y  c{,)  c  -  rU> 

i=l 


=  c%  CSAB  +  y  C2  CSAi  * 
i*l 


where  the  triple -sub script  coupling  coefficients  can  be  obtained  from 
the  variable  terms  in  Eq.  {17}  as 
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(51) 


CBAE‘  CWAB‘  CSAB  =  vB(b»  2  +  6)  +  vB(b’  y  -20o  +  8)‘ 

at  6  =  -  1,  *  +  eo,  -{Tr-eo}; 

CBAW»  CWAW’  ^SAW  =  vB^E»w'f'6o  "  ^  * 

at  0as  “  r* eo*  -(*-  e0)  j 

CBAS‘  CWAS’  CSAS  =  vB(pE*  *  "  6o  +  » 

at  0  =  ~  T  .  tt  +  9Q»  -  (tr  -  0O)  ; 

CBAi’  CWAi'  CSAi  =  VB(pt’I  +i0o+0)  +  vB(Pi'  "T  -(i+2)G0  +  e)  , 

at  e  =  -  —  (  w+  0O,  -(jt  -  ©0)  ; 

CAiB  =  vB(b«  f  -2i60-Gj  +  Vgjb,  ~  -  2(i+l)eo  -  j  . 

at  6  =  (j  -  i6o)  5 

CAiW  35  vB<pE’  *  +  {2i+i  )0o  +  6  )  .  .te-  f  -iG0; 

CAiS  ~  VB^pE’  W  -(2i+1)e0"S)  '  at  0=  -i0o  ; 

CAik  =  vB(pk,f-ieo-6)  +  vB(pk,^  -(3i+2)0o-eJ. 

at  8  =  I  -  i0o  5 

CAWA  =  VB^PE'  ®-6o»  UW}  at  9  =  0o  i 
CAWB  =  vB(Pe*  2lt  -6~  9o’  nW*  at  0  =  9o  ’ 
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(cont)  °ASA  =  VB^E*  °o  + 
and 


CASB  =  VB^PE»  6o  “ 


at  e  *  eQ; 

5,  at  e  *  . 


Equations  (49),  (50),  and  (51)  car  then  be  to  solve  any  order  of 
coupling  coefficients  from  £q.  (48). 


APPENDIX  D 

THE  SCATRAN  SOURCE  PROGRAMS  FOR 
NUMERICAL  COMPUTATION  BY  IBM  7094 


The  sourct  programs  fox  Eq.  (29)  and  the  associated  subroutine 
VB  presented  here  for  reference.  The  input  data  required  are  the 
horn  length  pE  in  terms  of  wavelength  and  the  angle  0j£  in  degrees. 

It  is  noted  that  the  notations  Q,  R,  rE»  and  cxE  in  the  programs 
correspond  to  A,B,  pE,  and  0£,  respectively,  in  Fig.  1. 


4? 


***** 

HOI 


H02 


H02A 


M02B 


integers  <h>- 

FLOATING  I NS  *  Nrf )  — 

COMPLEX  (CQRG.CGRI .SCURI .C»Q2»C*Ol <SC«Qi  I— 

COMPLEX  ( YR.XRS8. T3P  >  — 

COMPLEX  <C*Q.CQR.CO.XGS.X(iR«XO.xRS»»:fiQ.XR«XRSP»- 
COMPLEX  lXROP.XRP.XWQ.X*R»XQS8»  YQS*  YUR.VQ.uft.UT  )- 
COMPLEX  I  COM2*  Tl »T2»T3>“ 

COMPLEX  <CEXPl«>- 

OS MENS  ION  <R(200 )  *C3 <200 » .XQ (263 )  *XR<200> «XRP<2CQ> >- 
DIMENSION  (XQSB ( 200 ) vYG(200) )— 

L I  TER ALS  ( P I . 3 * l A  I  59265 . R AD . . 0 1 7*5329  *P 1 2 . 6* 2831053)- 
DIMENSION  (YR(200)tXRSfctC200) >- 
LITERALS  <COMZ«O..I.O. ;PIOH. | *57079631- 
1NPUT 

READ  INPUT  ♦8.INCSI- 
00  THROUGH  (H23 j . IN=0» l * |N*L*NCS- 
READ  INPUT  « HFO I  *  I RE »  ALPHE ) - 
CALCULATION  OF  CONSTANTS 

alr*rao*alphe- 

H*PIxI2»*AlR>- 

PROVIOEO  «H*E«PI/'<2*AALR»»tHxH-l- 
NS* <  2«*ALR  >/P|  — 

NW  =  1 2  •  -NS  )  - 
H=2*#RE#S|N«  <ALR>- 
8*R— 

do  through  <ho2>. i *1 « i * i •le.h- 
R{ I |xR( 1-1 )*C0S. (AlR)+R*COS. I I»AlR>- 
TRANSFER  (H02>  PROVIDED  (!.GE»H>- 

CALL  SUBROUTINE  IRVB.UVd » *VB* (PJ2.RE.P I0«- C I +1 >*AlR»2» »- 
C0<  l l*RVB+. I .UVB- 
CONTINUE  - 

CALL  SUBROUTINE  IRVE *UVBI«VB* IPJ2«6*PI-2.*AlR»2« >- 
C3RQcRv6+« I .UVB- 

CALL  SUBROUTINE  t Rve.UVB )  =  VB. (PI 2 .8.0. .2. ) - 
CGRQ«CORQ+ ( RVB+. I »UVB i— 

CALL  SUBROUTINE  (RVB.UVB » *VB. < PI2.RE.P1 0M-AlR«2. >- 
COR  1  iRVS+.  I  .UV8-* 

SCOR I *COMZ— 

00  THROUGH  Ih02A)*1*H-I*-1«I.GE*1- 

CALL  SUBROUTINE  (RVSI .UVBt >*VB» (Pi2*R( I >. I«ALR*2* >- 
CALL  SUBROUTINE  (RVB.UVB JcVU. (PI 2»R(  l I «PI - I  I *Z > »^lR«2» >- 
RVB I *RvB l ♦SVB- 
UV81 »UVBl*UVB- 

CALL  SUBROUTINE  < RVB.UVB > *VB. <P I  2 .RE.P I0H- < I f I J *ALR.2. ) - 
SCQR1  =SCQR  I  +  (RVBi  *  I  .UVB1  >*  (RVS+  •  I  .UVB }- 
CQR= I  COR I +SCOR I »✓<< !•♦• 1.0. I-CQRO)- 
CALL  SUBROUTINE  (RvB.UVB ) a VB. ( PI 2 »B. 3. »Pl OH+ALR.2, J- 
C*Q2»RVB+. I *UVS- 

CALL  SUBROUTINE  IRVB«UVB>*VB.CPI2.O.5.*PI0h-AlR»2»)- 
cwQ2«<c«Q2+cRva+*i*uva>  >*cqr- 
CALL  SUBROUTINE  (RvB.UVB )=vd. 1 P I 2. RE.P 12. 2. )- 
C*QI *RVB*« I .uva- 
SCWQt*C0MZ~ 

DO  THROUGH  (H02B ) . I sM1 1 »— l . I .»E« I - 

CALL  SUBROUTINE  (RV8 1 .UVB l )« VB. (PI 2.R ' J >. 3* «P10H* i ! 4 1 >*ALR 
.2*  )  — 

CALL  SUBROUTINE  (RVB.UVB I  * VB* < PI 2.R( I } .5. »P IOH- ( !♦» >*ALR«2 

• 

RVB1=RVB1-»RV8- 
UV8 I *UVB I 4UVB- 

CALL  SUBROUTINE  (RVB.UVB ) * VB. (PI 2.RE «P{ Oh- ( I ♦! } *AlR*2. I  — 
scwoi  *sc»oi  +  (Rvb+.  I  «uva>«{Rvai-* .  i  .uvbj  >- 
CW(3aC*Ql  ♦CW32+SCW0 1  — 

CALCULATION  OF  VARIABLES 


DO  THROUGH  (H22  )  <K>Di  J  «K.t_E»  1BQ- 
TH£TA=K*RAD- 

CALL  SU8R0UT I (RVB.UVB ) * VB. <PI 2.RE.P 1 -ALR+THETA.2. )- 
XQS'RVB-.*  I.UVB- 

CALL  SUBROUTINE  < RVB .UVB ) * VB. I P I 2 «R.P I0H+THETA.2 • >- 
XQR*RVB+. I .UVB- 

CALL  SUBROUT  5NE  (RVB.UVB ) * VB* IP  1 2  vR. 3. *Pl 0H-2«*ALR+TH£TA«2 

• 

XQR»XQR'MRVB+.I.UV6>- 

yQS*CEXPl* <  0  »  +  » U-P J2*RE*CQS. (THETA~AlR) >- 
YQfi*CEXPL.<0.-*-.  I.-PI2#R*SJN.  (THETA)  >- 
00  THROUGH  (H03 ) .  1 sh— i 4—1*1 *G£» 1  — 

CALL  SUBROUTINE  ( RVB *UVB > = VB. (P 12 «R< I ) .PI 0H+ l AALR+THET A.2» 

>- 

XQ(  2  >*RV6+.  t  .UVU- 

CALL  SUBROUTINE  < RVB.UVB ) =V3. ( P 12 . R ( I ) « 3. «P IOH- < I +2 >#ALR+T 
H£TA.2. )- 

H03  XQ< I ) =XQ< I >♦  (RVB*  »  I  .UVB  >- 

S=  !  •  - 

PROVIDED  (THETA. G.AlR) «S=0.- 

TRANSFER  (H05)  PROVIDED  (THETA. G.PIOH)- 

CALL  SUBROUTINE  (  RVB  .UVB  )  *V6.  (  P12  .RE.  PI -ALR-THETA.  2.  )- 
XRS»RVB+* 1 .UVB- 

CALL  SUBROUTINE  (RVB.UVB >=VB. (PI2.R.PI0H-THETA.2. )- 
XRQ*RVB+* I .UVB- 

CALL  SUBROUTINE  (RVB.UVB )=VB* (PI2«R.3**P10h-2.#AlR-THETA«2 
.  )  — 

XRQ=XRQ+(RVB+. l.UVB)- 
DO  THROUGH  (MO* > . I »H-» .-1 . 1 .GE. 1- 

CALL  SUBROUTINE  (RVB.UVB ) =VB» (PI2*R( I ) .PIOH4 I *ALR— THETA»2o 
l- 

XR( I )»RVB+. 1 .UVB— 

CALL  SUBROUTINE  (RVB  »UVB  >*VB.  (PI2»R(  I  )»3.  API  OH-  (  1 4-2  )*ALR— T 
META. 2. )- 

HO*  XR( 1 )*XR( 1 )+(RvB+. 1 .UVB>- 

TRANSFER  (H07>- 
H05  XRS3COM2— 

XRQ»C0M2- 

DO  THROUGH  <H06) » I =H-1 ,-I . I .GE. I- 
H06  XR< I }*COMZ— 

HOT  TRANSFER  (HOB)  PROVIDED  ( THETA. L.ALR )- 

CALL  SUBROUTINE  (RVB.UVB )  =  VB. (PI  2  »RE»  THET  A-AlS.N* )  — 
XwQxRVB-..  l.UVB- 

CALL  SUBROUTINE  < Rve «UVB ) *V8. (P I2.RE.P I 2-THETA-ALR.NW > - 
x«R=Rve+. i .uve- 
TRANSFER  (M09>- 
HOB  XWO«COMZ- 

XWR«C0M2~ 

M09  TRANSFER  (HI l )  PROVIDED  CK.L. 180.-ALPHE )- 

CALL  SUBROUTINE  5 RVB .UVB > «VB. <P I 2 .RE. 3.*pI -AlR-THETA.2* )- 
XRSP«RVB*. I *UVB- 

CALL  SUBROUTINE  (RVB.UVB ) « VB. (PI 2  *R«5.*P JOH-THETA.2. )— 
XROP*RVB+» I .  l  IVB— 

CALL  SUBROUTINE  (RVB .UVB } * VB. (PI2.R. 7. *PI0h-2,AAlR-THETA.2 
.  )- 

XRQP«XRQP+ («ve+. I .UVB  >- 
DO  THROUGH  (HlO).I*H— I»— I.I *GE* I— 

CALL  SUBROUTINE  ( RVB.UVB >■ VB. (P 12. R( I ) .S*AP!OH+ I aAlR-ThETA 

.2*  )~ 

XRPI I >.RVB+. I.UVB- 

CALL  SUBROUTINE  ( RVB.UVB )*  VB. (PI2»R( I ) »7»*P! Oh- ( I *2 )*AlR—T 
HETA.2. >- 

HlO  XRP (  I  ) ■XRS>{  I  )♦  (RV6+.  I  .UVB  )  — 

TRANSFER  (H!3>- 
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HI  1 

H  1 2 
H  i  3 


M 1 4 

NtS 


H 16 
H}9 


HI  9* 

H190 

C  ***** 


M20 

H2i 


C  ***** 


H2JA 


H22 


*R5P«C0*Z- 

XROP*COMZ- 

DO  THROUGH  (hi 2  >  *  I *H— l •- 1 • C  «G£* 1 - 
XRP< I )«COMZ- 

do  through  <hisj«i i«h-i.-i,i j.ge.i* 

TRANSFER  <H|4J  PROVIDED  < <.L.90»- (  *  J  ♦  I  >*ALPHE*0R*K«G.90.- I 
I  #ALPHE  )- 

YG<  1 1  »**CEXPl*  (O.  +  .  I  .  ( -P I  2*R  (  I  t  >*SINi  (  !  1#ALR*THETA)  )  }- 
CALL  SUBROUTINE  <RVB*UVB  )  *  VB.  IP  l  £*RE«P  I-(  2*  II  *i  )*ALR-THET  * 
»2»  I** 

XQSB( I  1 >»RVB** I •UVB- 
TRANSFER  (HI  5 )  — 

YC( I i i«COMZ- 
XOSB< I 1 I-COHZ- 
CONTiNUE  - 

TRANSFER  <H16  >  PROVIDED  (K ,L« (2*H*i J *AlPH£-1B0» .OR.K.G* (H+ 
i >*ALPH£-90» }- 

YR(H)*CeXPL. (O.  +  . I »  <-PI2#R(Hl*SlN.< (M-fl >*AlR*TH£TA>  J >- 
CALL  SUBROUTINE  (RVB«UVB ) * VB. (PI 2<R£«P I-( 2*H+ J J#ALR+THETA« 
2*>- 

xrsb<m»«rv8*«  i.uvb- 

TRANSFER  (H19S- 
XRSB(H)*C0M2> 

YR(Hl*COMZ- 
CONTINUE  - 

TRANSFER  (HI9AJ  PROVIDED  (<.L«90*~(H*l >*ALPME.OR.K,G6 J80»- 
(2*H*1 )*ALPHEI- 

CALL  SUBROUTINE  ( RVB.UV3 ) *VB. 5P 12 »RE. P l-( 2*H*J >*AlR-TH£TA* 
2*  >- 

XQSS  (  H  )  an  VB*  •  I  .  UVb- 

YCilHlsCEXPL. (0.*«1.(-P12*R(H)*SIN.( (H*l >*AL«-THETA >  >  }- 
TRANSFER  (Ml9t>>- 
XQSStHJaCOMZ- 
YG(H J»COMZ— 

CONTINUE  - 

CALCULATION  OF  OUTPUT  VALUES 
UR-S* YQS+XOS+  <  XRS+XRSP ) *YOR— 

T2P*C0NZ— 

Tl-COMZ- 

T2*COMZ— 

T3*CONZ- 

OO  THROUGH  (H20)t l«H-it-j«I*GE*l- 
T I *T l +CO( I 1*XG( I )- 
T2»T2+CQ( I )* (XR( l >*XRP( I  I  >- 
OO  THROUuf.  (H21  >  «  J  *H*  —  I  «  J  «GE« {— 

T3P«T3P*XGSB( 5  )*rO(  I  )- 
T3-T3+XRS B< I >*YR( l )- 

UT«UR*CWO*(XWQ*X*/R > *YGS*CQR*XGR*T I ♦( COR#  (XRQ4-XRQP  >+T2 }*TQ«+ 
T3+T3P- 

PROVIDED  (THETA, NE. 0. ) .  TPAnSFEP  (H2|A>~ 

URZ»CFA@S*CUR1- 

UTZ«CFABS*<UT?- 

PRJ NT-OUT  OF  VARIABLES 

I ALPHE * I O • • ALPHE - 

!RE*IO**«E- 

PUNCH  CARDS  »HFG6» < I ALPHE* IRE >- 
CONTINUE  - 

FuT»20»*FlOGI « {•ABS.OT/uTZj- 
FUR«20, *FLOG l • ( •  ABSeUR/'URZ)- 
IFUTIK  >«(FUT*90*  >*20»~ 

IFuR(X  >*<FUR*90. J*20»- 
KK (K } “ 6 • 666666 TAX- 

DIMENSION  ( IFUT( 181 )*  3FU«( IBS  >.**<181 » }- 
CONTINUE  - 
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H23  PUNCH  CAROS  *HFQ? * < <KK U > ♦ IFUT < I > ♦ I *0« I » I »LE» 190 >>- 

F  HFOl  < 2F >0*5) 

F  HFOS  <  2 (2X*  1  3  >  }  - 

F  HF07  <  S8I4 )- 

CAU.  SUBROUTINE  ( ) «ENDJOe. ( )- 
END  PROGRAM  IHOti- 


#*# 


SCATRAN 

SUBROUTINE  (RVd.UVB ) *VBo t r< «  ARG i ANGLE -FN 

provideo(  »Ads» ( .abs. angle-3.  i4isv265,'»'  e 

< ALTER*  - 

CONST  =  3. 14159265/FN  - 
COSANG  =  COS*  <  ANGLE )  - 
X  =  F<#ARG#i 1 ,+COSANG >  - 

call  function  <c> =sfresn. <x i  - 
CALL  FUNCTION  (  S  )  *UFRESN.  <  X  >  - 
RCOEFF  =  COS. <FK*ARG*COSAN-  J  - 
UCOEFF  =  SIN. (FK*ARG*COSANG ?  - 
RREST  -  .5-C/2.-S/2.  - 
UREST  *  S/Z.-C/2.  - 


RUVBOOOl 

)  - 

» .0001 ) .TRANSFER  TO 


*  *  */Fn>*SIN.  (CONST  )  *2.*.ABS.s  COS.  I  ANGLE/2.  )  >  /  t  COS.  (  CO 
NSTI-COS.  <  ANGLE/FN)  )  - 
RVB=COEF*(RCOEFF*RREST-UCOEFr*UREST  )  - 
UVB=COEF*<RCOEFF*UR£ST+UCOEFF*RREST )  - 
NORMAL  EXIT  _ 

ALTER  RVB*-.5«C0S. (FK#ARG )- 


FRES0001 


UV6=+.5*SIN. «F<#ARG)_ 

NORMAL  EXIT  - 
END  SUBPROGRAM  - 
END  PROGRAM  - 
***  STATEMENT  LISTING 
***  SCATRAn 

function  (o  =  rfresn.  <xs  - 

PROVIOEO  ( X. G* I  0. >  .  TRANSFER  To  (AsYmP)  - 
H  =  SORT. (2.#X  -3. 141S92?)  - 
C  a  H  - 

PROVIDED  I . ABS. CH ) .LB. .0000] ) •  TRANSFER  TO  (DONE)  - 
DO  THROUGH  (LOOP).  L= I «1. PROVIDED  ( .ABS. (H ) ,G. .00001 )  - 
H*w*< ( <3.-4.*L )#X*X )/( (4.*L+1 . >*(2.*L)«<2.*L-1 . > ) )  - 
C  *  C  T  H  »• 

LOOP  CONTINUE  - 

TRANSFER  TO  (DONE)  - 
ASYMP  si  =  I.  - 

AS1  =  SI  - 

STEMPI  =  .ABS. (SI)  - 

DO  THROUGH  (LOOP!).  L*2 ,2 .PROV I  DEO  ( .  Ao5« (SI )*G< .000 i * AND. ST 
EMPl / • ABS» { S| ) *6E* 1 • )  — 

STEMPI  s  o  ABS**  ( S  l  -5  - 

51  *  Sl*< <-| «* (2. *L-) • )*(2.*l-3. )  )XJ4»*x*X’ )  - 
AS  1  =  ASl  ♦  SI  - 

LOOPl  CONTINUE  - 

52  *  *  5/x  _ 

AS2  =  S2  - 

STEMP2  =  .AB$.(S2)  - 

DO  THROUGH  (UOOP2).  L'O  ..PROVIDED  C .ABS. (S2)«G*. 000] .AND. ST 
EMP2/.ABS.(S2),GE.l. )  - 
STEMP2  *  . ABS« (32  )  - 

S2  *  S2*t (-I «* (2«*L-3. )♦ (2.*L-1 • ) )/ (4. *X*X ) )  - 
AS2  =  AS2  ♦  S2  - 
L00P2  CONTINUE  - 

C  =  .5  +  <SIN»  (XI/SURT.  <6.283i834*X ;  )*ASl  -s-  (COS.  (X  )/SUR7 .  (  a 
»2S3I854*X) )»AS2  - 
DONE  NORMAL  EXIT  - 

END  SUBPROGRAM  - 
END  PROGRAM  - 
STATEMENT  LISTING 


#■#* 

«*t 


SCATRAN 

FUNCTION  !S)  =  UFRESN.  (X)  - 
PROVIDED  (X.S.lO. ).  TRANSFER  TO  (ASYMP)  - 


ufreoooi 
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LOOP 

ASYHP 

LOOPS 

LOOP2 

DONE 


H  =  SORT*  (2.*x/3*  14  Sb<i27  >*x/3»  - 

S  =  H  _ 

DO  THROUGH  < LOOP ) ♦  L  =  1 ♦ J *PROV ! DEO  ( . AdS. ( H > .G. .0000 J )  - 
H  =  H*  {  {  (  J  •— 4  •  *L  )*X*X)/(  (  <i »  *L+3»  )  *  { 2  »  *L+  1  . ) *  (2 »  *L  )  )  )  - 
S  s  S  H  - 
CONTINUE  - 

TRANSFER  TO  (DONE)  - 
SI  *  1.  - 
AS 1  *  SI  - 

STEKPl  •-  .A0S.IS1)  - 

00  THROUGH  {LOOPS).  L=2. 2. PROVIDED  ( . AdS. <S 1 > .G. . 000 J . AND.ST 
EMPJ/.ABS. <Sl S.GE.l . )  - 
STEMPl  =  .ABS. <SJ  )  - 

51  s  Sl*<  {-1  •  *  <2.*L-1  .  )*<2.*u-3.  )  )2(4.*X*X>  )  - 
ASt  *  ASS  +  SI  - 

CONTINUE  - 

52  s  .5/X  - 
AS2  *  S2  - 

STEHP2  *  .ABS*  t  S2 )  - 

OO  THROUGH  (LOOP2).  L*3« 2 .PROV I DED  { . ABS. (S2 > .G..OOOI . AND. ST 
£kPH/« ABS. I S2 ) »GE . I « )  — 

STEMP2  -  .A0S.IS2)  - 

S2  *  S2*< (-I .<M2.*L-3. )*C2.»L-I • ) >/{4.*X*X) )  - 
AS2  =  ASH  +  52  - 
CONTINUE  - 

S  *  .5  -  (COS. (XJ/SCIRT «  (6»2831854*X ) ) *AS1  -  (SIN. (Xj/SOST. (& 

•283ia54*X) )*AS2  - 

NORMAL  EXIT  - 

END  SUBPROGRAM  - 

END  PROGRAM  - 
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